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Abstract 



Three dimensional implementations of liquid state theories offer an efficient 



alternative to computer simulations for the atomic-level description of aqueous 
solutions in complex environments. In this context, we present a (classical) 
O molecular density functional theory (MDFT) of water that is derived from first 

principles and is based on two classical density fields, a scalar one, the particle 

43 

density, and a vectorial one, the multipolar polarization density. Its implementa- 
tion requires as input the partial charge distribution of a water molecule and three 
measurable bulk properties, namely the structure factor and the k-dependent lon- 

00 

gitudinal and transverse dielectric constants. It has to be complemented by a 

solute-solvent three-body term that reinforces tetrahedral order at short range. 

The approach is shown to provide the correct three-dimensional microscopic sol- 
t— ( vation profile around various molecular solutes, possibly possessing H-bonding 

sites, at a computer cost two-three orders of magnitude lower than with explicit 
^ simulations. 
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The numerical methods that have emerged in the second part of the last century 
from liquid-state theories (lj|2), including integral equation theory in the interaction- 
site [3- 8 1 or molecular (9 14 picture, classical density functional theory (DFT) [15 17 



or classical fields theory 18 21 , have become methods of choice for many physical 



chemistry or chemical engineering applications 22 25 . They can yield reliable predic 



tions for both the microscopic structure and the thermodynamic properties of molecular 
fluids in bulk, interfacial, or confined conditions at a much more modest computational 
cost than molecular-dynamics or Monte-Carlo simulations. A current challenge con- 
cerns their implementation in three dimensions in order to describe molecular liquids, 
solutions, and mixtures in complex environments such as atomistically-resolved solid 
interfaces or biomolecular media. There have been a number of recent efforts in that 



direction using 3D-RISM 26 31 , lattice field 32, 33 or Gaussian field 20,34 theo 



ries. Recently, a molecular density functional theory (MDFT) approach to solvation 



has been introduced. 35 40 It relies on the definition of a free-energy functional de- 



pending on the full six-dimensional position and orientation solvent density. In the 
so-called homogeneous reference fluid (HRF) approximation, the (unknown) excess 
free energy can be inferred from the angular- dependent direct correlation function of 
the bulk solvent, that can be predetermined from molecular simulations of the pure 
solvent. Compared to reference molecular dynamics calculations, such approximation 



was shown to be accurate for polar, non- hydrogen bonded fluids 35,38-41 , but to 



require some corrections for water 39 42 44 



In this paper we introduce a simplified version of MDFT that can be derived 
rigorously for SPC- or TIPnP-like representations of water, involving a single Lennard- 
Jones interaction site and distributed partial charges. In that case we will seek a simpler 
functional form expressed in terms of the particle density n(r) and site-distributed 
polarisation density P(r), and requiring as input simpler physical quantities than the 
full position and orientation-dependent direct correlation function. This functional 
may be considered as a multipolar generalization of the generic dipolar fluid free-energy 



functional J-[n(r), P(r)] introduced in Refs. 35,37 



A functional o/n(r) and P(r). We start from SPC- or TIPnP-like representation 



of water, constituted by a single Lennard- Jones center, located on the oxygen, and m 
charges distributed on various sites. Each molecule is supposed to be rigid with position 
r-j and orientation For a given water configuration, we define the microscopic 
particle densities 

N 

p(r,n) = ^25(T-Ti)5(n-ni), (l) 



i=l 
N 



Pn{r) 

1=1 



- ^5(r- ri )= f dQp(r,Q), (2) 

i=l J 

and the charge and multipolar polarization density 

N 

p c (r) = y]<r(T — T h Sli) 

i=l 

= J dr'dQa{r-r',tt)p(r',tt), (3) 

JV 

P c (r) = ^ M (r-r,,a) 

i=i 

= fMHr-S.aW,*). (4) 



cx(r, Q) is the molecular charge density and, according to the definition of Refs 46,48 , 
fi(r, Q) is the molecular polarization density of a water molecule taken at the origin 
with orientation Q: 

a(r,Q) = J^g OT 5(r-s m (fi)), (5) 

m 

fj,{r,Q) = y2q m s m (Q) / du8(r- us m (tt)), (6) 

J ° 

where s m (fi) indicates the location of the m th atomic site for a given Q. It can be 
easily checked that molecular charge and polarization densities are linked by the usual 
relation <r(r, Q) = — V ■ n(r, Q). In k-space 

M (k,n) = -i^^M-^^-i) (7) 
= ij,(n) + z -J2<im(k-s m (n))s m (n) + ..., (8) 

m 
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with = J^mlrnSj^Q) — fj, ft (/io the molecule dipole; for water, f2 is the unit 

vector along the O-H bonds angle bissector), such that the molecular polarization 
density reduces to a molecular dipole located at the origin at dominant order, but does 
include the complete multipole series at all other orders. 

Our derivation starts from the observation that the Hamiltonian of N water 
molecules in the presence of an embedded solute, described by an external molecular 
force field, can be written as 



where T and U are the water kinetic and pair-wise potential energy, respectively. 
$ n (r), E c (r) denote the value of the external Lennard- Jones potential and electric 
field at position r 



The solute is here described by atomic sites j, located at Rj, with Lennard- Jones 
parameters (Jj,£j (using Lorentz-Berthelot mixing rules with respect to water), and 
point charges Qj. 

Following the original derivation of Evans [15] and extending it to four indepen- 
dent external field variables, $n(r) and E x ^ tZ (r), instead of just one, it can be easily 
proved that the grand potential of the solute-water system at a given water chem- 
ical potential may be expressed as a functional of the one-particle number density, 
n(r) = (p n (r)), and of the one-particle polarization density, P(r) = ^P c (r)^>, that is 
= 0[n(r), P(r)]. Minimization of this functional with respect to those two fields 
yields the equilibrium densities and the value of the grand potential. Taking as a ref- 
erence the bulk water system at the same chemical potential, with number density no 
and grand potential 0o, the same properties hold true for the solvation free energy 
J-"[n(r), P(r)] = 0[n(r),P(r)] — ©o- This constitutes the first important result of this 
work. 




(9) 




(10) 



(11) 
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Expression of the functional. We proceed by invoking an equivalent of the Kohn- 
Sham decomposition in electronic DFT. Introducing the intrinsic Helmholtz free energy 



JF[n(r),P(r)]=JF int [n(r),P(r)] + J dr $ n (r)n(r) - / dr 



E,W-P(r) 



:i2i 



and decomposing it into ideal (non-interacting) part and excess part, Ti nt = J^id + J^exc, 
we go back to the full one-particle angular density p(r, fl) = (p(r,Q)) to express the 
ideal term (equivalent to going back to orbitals in eDFT) 

'87r 2 p(r,fi)~ 



Fid[p] = k B T J drdVt 



p (r, Q) In 



no 



p(r,n)) + 



n 



57T 



(13) 



where dVt indicates an angular integration over the three Euler angles. The spatially 
and orientationally uniform fluid density no/87r 2 is used as reference. We also expand 
the excess term around the homogeneous liquid state at that density 



^[n(r),P(r)] 



= T int [n{r), P(r)] - F id [n{r), P(r)] 

= J dv 1 dv 2 S~ 1 (r 1 2) An(ri) An(r 2 ) 

+ -i- / dr 1 dr 2XL 1 (r 12 )P L (r 1 )-P i (r 2 ) 
ovre J 

+ ^— I rfr 1 rfr 2 x T 1 (ri 2 )PT(ri) • P T (r 2 ) 



(14) 



- k u T I dr An ^ - k B T I <lr 



■PW 



J- cor [n(r),P(r)]. 



2n 7 2/ign 

The first three terms stem from the quadratic expansion of Tint- Here Pl, Pt designate 
the longitudinal anrd transverse components of the polarization vector, that are defined 
in k-space by 

p L (k) = (p(k) • k)k, p T (k) = p(k) - p L (k). (is) 

5'~ 1 (r) and xZr( r ) are defined as the Fourier transforms of the inverse of the structure 
factor S(k) and longitudinal and transverse electric susceptibility Xi{k) and Xr{k), 
which may be obtained either from experimental data or by simulations through the 
correlation of the molecular density fluctuations and of the longitudinal and transverse 



polarisation vector fluctuations 45-48 . We have computed these quantities for SPC/E 
water [49] using the k-space method described by Bopp et al 47,48 ; see Fig. [l] The 
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dielectric susceptibilities are related to the longitudinal and transverse dielectric con- 
stants through Xi{k) = 1 — l/ez,(fc) and Xt(&) = ( e r(&) — l)/47r. Their very small k 
behavior, impaired in simulation by finite size effects, can thus be extrapolated using 



the macroscopic value of the dielectric constants ex(0) = 6t(0) = 71 50 



Note that in eqn 14 the longitudinal polarization term can be also written as a 
charge-charge interaction term involving the bound charge density p c (r) = — V • P(r) 
(that would also include the free charges in the case of a ionic solution). Note also 
that we have neglected the density-polarisation correlations which we verified to be 
very small in water. The negative terms in eqn [14] stand for the linearized translational 
and orientational entropy, whereas the last term is the unknown correction term that 
contains supposedly all orders in An(r) and P(r) higher that 2. 

As will be seen below, the functional with J-" cor = tends to underestimate the 
tetrahedral order in the vicinity of the solute atomic sites giving rise to H-bonds with 
the solvent. Such defect was also present in our previous MDFT functional, involving 



the full angular-dependent c-functions in the excess functional 39 . To correct it 



we introduce as before a short-range three-body potential term borrowed from the 



coarse-grained water model of Molinero et al. for water and ions in water 51 , 52 that 
enforces some tetrahedral order around all atomic sites k susceptible to give a H-bond 
interaction 

Fcar[n(r)] = ^fceT^tAfc / dr 2 dr 3 f k (r k2 ) /fc(rjfc 3 )( — — - cos6> ) 2 n(r 2 ) n(r 3 ) 

£ , J T k2 T k3 

k 

+ \\ I dr 2 dr 3 f k (r k2 ) fw^) ^ 2 — - cosflo) 2 n(r 2 ) n(r 3 )], (16) 

J r fc2^23 

with r ki = |r.j — H k \. The two terms describe three-body H-bonding to the the first and 
second water solvation shell, respectively. f k (r) and f w (r) are the short-range functions 
introduced by DeMille and Molinero for ion-water and water-water contributions. Their 
systematic parametrization, as well as the choice of the strength parameters A^., A 2 ,, 
were alluded to in Ref. [39]. For simplicity, we have selected below the values /c^TA^. = 
fcfiTA^ = 75 kJ/mol for ionized H-bonding sites, and 50 kJ/mol for neutral ones, 
whatever their chemical nature. 
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We note finally that the above functional 12 14 (with J-" cor = 0) can be shown 
to be equivalent to the generic dipolar functional introduced in Refs [35] [37] if the 
molecular polarization density is restricted to dipolar order in eqn[8j such that P(r) = 
p J dQflp (r, Q). The present work can thus be viewed as a generalization of Refs 
37] to a Stockmayer-like multipolar solvent. 
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Implementation and results. The functional defined by eqs 12 to 14 can be mini- 
mized in the presence of the external potential by discretizing p (r, Q) on a position and 
orientation grid. The minimization is performed in fact with respect to a "fictitious 
wave function", \I/(r, Q) = p(r,f2) 1//2 , in order to prevent the density from becoming 
negative in the logarithm. The position are represented on a 3D-grid with typically 
3-4 points per Angstrom whereas the orientations are discretized using a Lebedev grid 
(yielding the most efficient angular quadrature for the unit sphere [53]) for 9, (ft, the 
orientation of the molecule C2v axis, and a regular grid for the remaining angle ip (from 
to 7r). We use typically Nq = 12 x 3 angles. The calculation begins with the tabu- 
lation of the external Lennard- Jones potential n ( r ), and external electric field E c (r). 
The latter is computed by extrapolating the solute charges on the grid and solving 
the Poisson equation. The molecular dipole function ^t(k, Q) is also precomputed in k- 
space. At each minimization cycle, p (r, Q) is Fourier transformed to get p(k, Q), which 
is integrated over angles to get n(k) = J dQ p(k, Q) and P(k) = f dQ fi(k, Q)p(k, Q), 
as well as its longitudinal and transverse components defined by eqn [15} The non-local 



excess terms in 14 are computed in k-space and so are the gradients with respect to 
p(k, Q); those are then back-transformed to r-space. The three-body contribution can 
be expressed as a sum of one-body terms and computed efficiently in r-space. For 



minimization, we used the L-BFGS quasi-Newton optimization routine 54 which re- 
quires as input at each cycle free energy and gradients. The calculation are performed 
on a standard workstation or laptop. The convergence turns out to be fast, requiring 
at most 25-30 iterations, so that despite the important number of FFT's be handled 
at each step (2 x Nq) [55], one complete minimization takes only a few minutes on 



a single processor. The computation of identical quantities using molecular dynamics 



simulations requires 10 — 10 times more computer time 40,41 . 
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To describe water, we have chosen the SPC/E model (49] which defines both the 



external potential in eqs [TOpl] and the structure factor and dielectric susceptibilities 
of equation [14J Functional minimization was performed in the presence of different 
solutes placed the center of cubic box of size L = 25 A, with 80 3 grid points. In 
figure 2, we display the ion-oxygen pair distribution functions obtained by MDFT 
minimization around two representative monovalent ions and compare them to the 
corresponding curves generated by molecular dynamics simulations using a similar box 
size. In the straight homogeneous reference fluid approximation, F corr = in eqn 14 , it 



appears as in Ref. 39 that the first peak exhibits the correct position and width but 
is overestimated in height, whereas the second peak is misplaced; this later deficiency 
points to a lack of tetrahedral order. Those features are well corrected by inclusion 



of the three-body term of eqn 16 Identical comparisons are carried out in figure 3 
for a SPC/E water molecule dissolved in our MDFT water model, and in figure 4 for 
a N-methyl-acetamide molecule (NMA: CH3NHCOCH3), the paradigm for a peptide 
bond chemical motif. Three-body contributions were added for the O and N sites with 
the parameters given in the previous section, the H-atoms being included implicitly 
in the associated heavy atom. The MDFT approach is able to provide not only pair 
distribution functions, but also the full three-dimensional solvent structure; this is 
illustrated for each molecule by an image of the surrounding three-dimensional water 
density. For water in water, the expected tetrahedral, four-fold coordination appears 
clearly, as does the fine structuring of water around the C=0 and N-H bonds of NMA. 

The purpose of this note was to derive a refined version of MDFT that is based 
on two natural physical fields, the particle density and polarization density, and on 
three simple, measurable bulk water properties. It was shown also that an efficient 
three-dimensional numerical implementation can be carried out. The next step will be 
the application of this new formalism to the hydration of complex molecular solutes, 
such as the atomistically-resolved clay surfaces considered in Ref. [4l] with a reduced 
dipolar water model, or to the hydration of biological molecules. We will focus also 
on a systematic study of the performance of our novel MDFT scheme for the accurate 
prediction of molecules hydration free energies, in addition to the microscopic hydration 
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structure. 

Acknowledgments 

The authors acknowledge financial support from the Agence Nationale de la Recherche 
under grant ANR-09-SYSC-012. DB acknowledges many thorough discussions about 
DFT with Rosa Ramirez and Shuangliang Zhao. 

References 

[1] Hansen, J. P.; McDonald, I. R. Theory of Simple Liquids; Academic Press: Lon- 
don, 1989. 

[2] Gray, C. G.; Gubbins, K. E. Theory of Molecular Fluids, Volume 1: Fundamentals; 
Clarendon Press: Oxford, 1984. 

[3] Chandler, D.; Hendersen, H. Optimized Cluster Expansions for Classical Fluids 
II- Theory of Molecular Liquids. J. Chem. Phys. 1972, 57, 3275-3276. 

[4] Hirata, F.; Rossky, P. J. An Extended RISM Equation for Molecular Polar Fluids. 
Chem. Phys. Lett. 1981, 83, 329-334. 

[5] Hirata, F.; Pettitt, B. M.; Rossky, P. J. Application of an Extended RISM Equa- 
tion to Dipolar and Quadrupolar Fluids. J. Chem. Phys. 1982, 77, 509-520. 

[6] Reddy, G.; Lawrence, C. P.; Skinner, J. L.; Yethiraj, A. Liquid State Theories for 
the Structure of Water. J. Chem. Phys. 2003, 119, 13012-13016. 

[7] Dyer, K. M.; Perkyns, J. S.; Pettitt, B. M. A Site-Renormalized Molecular Fluid 
Theory. J. Chem. Phys. 2007, 127, 194506-194506-14. 

[8] Dyer, K. M.; Perkyns, J. S.; Stell, G.; Pettitt, B. M. A Molecular Site-Site Integral 
Equation that Yields the Dielectric Constant. J. Chem. Phys. 2008, 129, 104512- 
104512-9. 



8 



[9] Blum, L.; Torruella, A. J. Invariant Expansion for 2-Body Correlations I- Thermo- 
dynamic Functions, Scattering, and Ornstein-Zernike Equation. J. Chem. Phys. 
1972, 56, 303-310. 

[10] Blum, L. Invariant Expansion II- Ornstein-Zernike Equation for Nonspherical 
Molecules and an Extended Solution to Mean Spherical Model. J. Chem. Phys. 
1972, 57, 1862-1869. 

[11] Patey, G. N. Integral-Equation Theory for Dense Dipolar Hard-Sphere Fluid. Mol. 
Phys. 1977, 34, 427-440. 

[12] Fries, P. H.; Patey, G. N. The Solution of the Hypernetted- Chain Approxima- 
tion for Fluids of Nonspherical Particles - a General Method with Application to 
Dipolar Hard-Spheres. J. Chem. Phys. 1985, 82, 429-440. 

[13] Richardi, J.; Fries, P. H.; Krienke, H. The Solvation of Ions in Acetonitrile and 
Acetone: A Molecular Ornstein-Zernike Study. J. Chem. Phys. 1998, 108, 4079- 
4089. 

[14] Richardi, J.; Millot, C; Fries, P. H. A Molecular Ornstein-Zernike Study of 
Popular Models for Water and Methanol. J. Chem. Phys. 1999, 110, 1138-1147. 

[15] Evans, R. Nature of the Liquid- Vapor Interface and Other Topics in the Statistical- 
Mechanics of Nonuniform, Classical Fluids. Adv. Phys. 1979, 28, 143-200. 

[16] Evans, R. In Henderson, D., Ed., Fundamental of Inhomogeneous Fluids, New 
York, 1992. Marcel Dekker. 

[17] Wu, J.; Li, Z. Density Functional Theory for Complex Fluids. Ann. Rev. Phys. 
Chem. 2007, 58, 85-112. 

[18] Chandler, D. Gaussian Field Model of Fluids with an Application to Polymeric 
Fluids. Phys. Rev. E 1993, 48, 2898-2905. 

[19] Lum, K.; Chandler, D.; Weeks, J. D. Hydrophobicity at Small and Large Length 
Scales. J. Phys. Chem. B 1999, 103, 4570-4577. 

9 



[20] ten Wolde, P. R.; Sun, S. X.; Chandler, D. Model of a Fluid at Small and Large 
Length Scales and the Hydrophobic Effect. Phys. Rev. £2002, 65, 011201-011201- 
9. 

[21] Coalson, R. D.; Duncan, A. Statistical Mechanics of a Multipolar Gas: A Lattice 
Field Theory Approach. J. Phys. Chem. 1996, 100, 2612-2620. 

[22] Gray, C. G.; Gubbins, K. E.; Joslin, C. J. Theory of Molecular Fluids, Volume2: 
Applications; Clarendon Press: Oxford, 2011. 

[23] Ravikovitchand, P.; Neimark, A. Density Functional Theory Model of Adsorption 
on Amorphous and Microporous Silica Materials. Langmuir 2006, 22, 11171— 
11179. 

[24] Gor, G. Y.; Neimark, A. Adsorption-Induced Deformation of Mesoporous Solids: 
Macroscopic Approach and Density Functional Theory Langmuir 2011, 27, 6926- 
6931. 

[25] Wu, J. Density Functional Theory for Chemical Engineering: From Capillarity to 
Soft Materials. AIChE Journal 2006, 52, 1169-1193. 

[26] Beglov, D.; Roux, B. An Integral Equation to Describe the Solvation of Polar 
Molecules in Liquid Water. J. Phys. Chem. B 1997, 101, 7821-7826. 

[27] Kovalenko, A.; Hirata, F. Three-Dimensional Density Profiles of Water in Contact 
with a Solute of Arbitrary Shape; a RISM Approach. Chem. Phys. Lett. 1998, 
290, 237-244. 

[28] F. Hirata, E. Molecular Theory of Solvation; Kluwer Academic Publishers, Dor- 
drecht, 2003. 

[29] Yoshida, N.; Imai, T.; Phongphanphanee, S.; Kovalenko, A.; Hirata., F. Molecular 
Recognition in Biomolecules Studied by Statistical-Mechanical Integral-Equation 
Theory of Liquids. J. Phys. Chem. B 2009, 113, 873-886. 



10 



[30] Kloss, T.; Kast, S. M. Treatment of Charged Solutes in Three-Dimensional Integral 
Equation Theory. J. Chem. Phys. 2008, 128, 134505-134505-7. 

[31] Kloss, T.; Heil, J.; Kast, S. M. Quantum Chemistry in Solution by Combining 3d- 
Integral Dquation Theory with a Cluster Embedding Approach. J. Phys. Chem. 
B 2008, 112, 4337-4343. 

[32] Azuara, C; Lindahl, E.; Koehl, P. Pdb .hydro: Incorporating Dipolar Solvents 
with Variable dDnsity in the Poisson-Boltzmann Treatment of Macromolecule 
Electrostatics. Nucleic Ac. Res. 2006, 34, W38-W42. 

[33] Azuara, C; Orland, H.; Bon, M.; Koehl, P.; Delarue, M. Incorporating Dipolar 
Solvents with Variable Density in Poisson-Boltzmann Electrostatics. Biophys. J. 

2008, 95, 5587-5605. 

[34] Varilly, P.; Patel, A. J.; Chandler, D. An Improved Coarse- Grained Model of 
Solvation and the Hydrophobic Effect. J. Chem. Phys. 2010, 134, 074109-074109- 
15. 

[35] Ramirez, R.; Gebauer, R.; Mareschal, M.; Borgis, D. Density Functional Theory 
of Solvation in a Polar Solvent: Extracting the Functional from Homogeneous 
Solvent Simulations. Phys. Rev. E 2002, 66, 031206-031206-8. 

[36] Ramirez, R.; Mareschal, M.; Borgis, D. Direct Correlation Functions and the 
Density Functional Theory of Polar Solvents. Chem. Phys. 2005, 319, 261-272. 

[37] Ramirez, R.; Borgis, D. Density Functional Theory of Solvation and its Relation 
to Implicit Solvent Models. J. Phys. Chem. B 2005, 109, 6754-6763. 

[38] Gendre, L.; Ramirez, R.; Borgis, D. Classical Density Functional Theory of Sol- 
vation in Molecular Solvents: Angular Grid Implementation. Chem. Phys. Lett. 

2009, 474, 366-370. 

[39] Zhao, S.; Ramirez, R.; Vuilleumier, R.; Borgis, D. Molecular Density Functional 
Theory of Solvation: From Polar Solvents to Water. J. Chem. Phys. 2011, 134, 
194102-194102-13. 

11 



[40] Borgis, D.; Gendre, D.; Ramirez, R. Molecular Density Functional Theory: Ap- 
plication to Solvation and Electron- Transfer Thermodynamics in Polar Solvents. 
J. Phys. Chem. B 2012, 116, 2504-2512. 

[41] Levesque, M.; Marry, V.; Rotenberg, B.; Jeanmairet, G.; Vuilleumier, R.; Borgis, 
D. Solvation of Complex Surfaces via Molecular Density Functional Theory. J. 
Chem. Phys. 2012, 137, 224107-224107-8. 

[42] Zhao, S.; Jin, Z.; Wu, J. A new theoretical method for rapid prediction of solvation 
free energy in water. J. Phys. Chem. B 2011, 115, 6971-6975. 

[43] Zhao, S.; Jin, Z.; Wu, J. Correction to New Theoretical Method for Rapid Predic- 
tion of Solvation Free Energy in Water. J. Phys. Chem. 5 2011, 115, 15445-15445. 

[44] Levesque, M.; Vuilleumier, R.; Borgis, D. Scalar Fundamental Measure Theory 
for Hard Spheres in Three Dimensions. Application to Hydrophobic Solvation. J. 
Chem. Phys. 2012, 137, 034115-034115-9. 

[45] Raineri, F. O.; Resat, H.; Friedman, H. L. Static Longitudinal Dielectric Function 
of Model Molecular Fluids. J. Chem. Phys. 1992, 96, 3068-3084. 

[46] Raineri, F. O.; Friedman, H. L. Static Transverse Dielectric Function of Model 
Molecular Fluids. J. Chem. Phys. 1993, 98, 8910-8918. 

[47] Bopp, P. A.; Kornyshev, A. A.; Sutmann, G. Static Nonlocal Dielectric Function 
of Liquid Water. Phys. Rev. Lett. 1996, 76, 1280-1283. 

[48] Bopp, P. A.; Kornyshev, A. A.; Sutmann, G. Frequency and Wave- Vector Depen- 
dent Dielectric Function of Water: Collective Modes and Relaxation Spectra. J. 
Chem. Phys. 1998, 109, 1939-1958. 

[49] Berendsen H. J. C; Grigera J. R.; Straatsma T. P. The Missing Term in Effective 
Pair Potential. J. Phys. Chem. 1987, 91, 6269-6271. 

[50] Kusalik P. G.; Svishchev I. M. The spatial Structure in Liquid Water Science 
1994, 265, 1219-1221. 

12 



[51] Molinero, V.; Moore, E. B. Water Modeled as an Intermediate Element between 
Carbon and Silicon. J. Phys. Chem. B 2009, 113, 4008-4016. 

[52] DeMille, R. C; Molinero, V. Coarse-Grained Ions without Charges: Reproducing 
the Solvation Structure of NaCl in Water. J. Chem. Phys. 2009, 131, 034107- 
034107-16. 

[53] Lebedev, V. I.; Laikov D. N. A Quadrature Formula for the Sphere of the 131st 
Algebraic Order of Accuracy. Dokl. Math. 1999, 59, 477-481. 

[54] Byrd, R. H.; Lu, P.; Nocedal, J. A Limited Memory Algorithm for Bound Con- 
strained Optimization. SIAM J. Scient. Stat. Comp. 1994, 16, 1190-1208. 

[55] Frigo M. and Johnson S. G. FFTW: a free collection of fast 
fourier transforms routines in C. Available at the FFTW home page: 
http : / /theory . lcs . mit . edu/ ~f f tw 



13 







10 

o -1 

k(A ) 



15 



20 



Figure 1: Structure factor and longitudinal and transverse dielectric susceptibility for 



SPC/E water computed according to the formulas in Ref 48 
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Figure 2: Ion-oxygen pair distribution functions for chloride and potassium ions in 
SPC/E water computed by MD (red lines) or MDFT without (blue lines) or with the 



three body term of eqn 16 (black lines) 
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Figure 3: Top: Isosurface of high water density (n(r) = 3rio) around a tagged SPC/E 
water molecule as determined by MDFT. Bottom: Corresponding site-oxygen pair 
distribution functions computed by MD (red lines) or MDFT (black lines). 
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Figure 4: Top: Three-dimensional water density around a N-methyl-acetamide 
molecule obtained by MDFT. Bottom: A selection of solute site-oxygen pair distri- 
bution functions computed by MD (red lines) or MDFT (black lines). 



17 



